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A covariant model of elastic pion-nucleon scattering based on the Bethe-Salpeter equation is presented. 
We obtain a good description of the S- and P-wave phase shifts up to 360 MeV laboratory energy. We 
also compare results from the if-matrix approach and several 3-dimensional quasipotential equations 
, to the Bethe-Salpeter equation. 

o 

p | 1 Introduction 

■ Over the past decade many dynamical models of pion-nucleon (irN) scattering based on meson- 

exchange have been developed. These models invariably begin with an effective hadronic Lagrangian 
describing the couplings between the various mesons and baryons. The tree-level diagrams obtained 
from this Lagrangian are then unitarized using an approximation to the Bethe-Salpeter equation 
(BSE) jjj, such as the iT-matrix approach or a 3-dimensional (3-D) quasipotential equation. The 
■^j- | 3-D reduction procedure is ambiguous because there are an infinite number of different quasipoten- 

tial equations [||, [|, each having different off-shell behaviour There is no overwhelming reason 
to choose one quasipotential equation over another. Sometimes quasipotential equations have been 
devised so as to have the correct one-body limit, however this argument has been shown to be not 
applicable to irN scattering ||. Furthermore, many of the commonly used quasipotential equa- 
tions violate charge-conjugation symmetry 0. Here we avoid these problems by constructing a 
covariant model of elastic ttN scattering Jjj] based on the BSE, which we solve without making any 
approximations to the relative-energy dependence of the kernel. 
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2 The model 

The BSE for the vriV -> ttN amplitude T is 



T(q',q;P)=V(q',q;P)-j^r i J d 4 q"V(q',q";P)G 7rN (q";P)T(q",q;P), (1) 

where G n N is the 2-body 7riV propagator. In principle, both the tt and iV propagators should be 
fully dressed, but for only 2-body unitarity to be maintained, G^n has the simple form 

G nN (q- P) = - n N P+j + m N 

(u, n P - q) 2 -ml + ie (^nP + q) 2 - m 2 N + ie ' 

where and ^ are functions of s = P 2 such that /U v + = 1. The solution of the BSE does 
not depend on the choice of (j-n(s) and /%(s), so we use the simplest possibility: u. N = fj, n = 1/2. 
The interaction kernel V is truncated to include only the 2nd order 7riV — > irN diagrams obtained 
from the following interaction Lagrangian: 

Ant = N 75 1 11 tN • cVr + C nNA + g uNN NNa+ 9 ^ a d^-Wir 

+9 p nnN \ t ■ (j^p " + -^L-o^&p ^ JV + g^p » ■ (tt x O^tv) . (3) 



Due to the ambiguity in the description of spin-3/2 particles, we consider two different possibilities 
for the 7riVA vertex: 



£na = -^ E ^e^(d lx A l/ ) l5la TN-dpir + h.c. (5) 

Here is the conventional ttNA vertex which contains the so-called off-mass-shell parameter 

xa, while £^y A is the Pascalutsa vertex [||. For the A propagator we use the standard Rarita- 
Schwinger (RS) form |J. It is well known that the RS spin-3/2 propagator contains off-mass-shell 
spin-1/2 components. The use of the Pascalutsa 7riVA vertex, together with the RS A propagator, 
ensures that the s- and u-channel A poles present in V are free of any contributions from the 
spin-1/2 components of the RS A propagator. 
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Table 1: The coupling constants and particle masses resulting from fits to the ttN data using the 
two different form factor parameterizations (denoted as model I and model II) , and the two different 
ttNA vertices. The quantities in boldface were varied in the fits. All masses are in GeV. 

Regularization is achieved by the introduction of form factors. We consider two different param- 
eterizations: in model I we associate a cutoff function with each vertex, where this cutoff function 
is taken as the product of form factors that depend on the 4-momentum squared of each particle 
present at the vertex. In model II, we associate a form factor only with the pion propagator. In 
both models, each form factor is chosen to have the form: 

This form factor fulfils the requirement of only having poles along the real axis. In this work we 
use ma = 1 for model I, and = 8 for model II. 



We solve the BSE by first expanding the nucleon propagator in the ttN intermediate states into 
positive and negative energy components, and then sandwiching the resulting equation between 
Dirac spinors. This gives two coupled 4-D integral equations which are reduced to 2-D integral 
equations after partial wave decomposition. A Wick rotation [|l0) is performed in order to obtain 
equations suitable for numerical solution. This means that all amplitudes are analytically continued 
in the relative-energy variables from the real axis to the imaginary axis, thereby avoiding the 
singularities of the kernel. Form factors with poles only along the real axis do not interfere with 
the Wick rotation provided the cutoff masses are large enough 0. 

3 Numerical results 

The free parameters are determined in \ 2 fits to the S- and P-wave phase shifts up to 360 MeV pion 
laboratory energy from the SM95 partial wave analysis Jll[]. The parameters obtained are shown 
in Table [[]. Note that g^NN was fixed at g 2 ^^/^ = 13.5, while g^N an< ^ m ff were determined 
by the nucleon renormalization procedure 0. This ensures that in the Pu partial wave the dressed 
s-channel nucleon pole diagram has a pole at the physical nucleon mass with a residue related to 
the physical irNN coupling constant. 

The resulting phase shifts are shown in Figure 1. When the conventional 7T-/VA vertex is used, 
we obtain very good agreement with the partial wave analysis. There is some disagreement in the 
Pu partial wave when the the Pascalutsa ttNA vertex is used, which suggests that the spin- 1/2 
components of the RS A propagator are necessary in order to obtain a good fit to the irN data. 
Note that the model II results are very close to the model I results, even though model II has four 
less free parameters. Therefore in our framework the use of a different cutoff mass for each particle 
results in unnecessary free parameters. The scattering lengths and volumes are shown in Table [||, 
where it is seen that we obtain reasonable agreement with the partial wave analyses. 

Our coupling constants are consistent with the commonly accepted values in the literature, 
with the exception of f-^NA when the Pascalutsa ttNA vertex is used. In this case, the irNA 
coupling constant is around twice as large as the so-called "empirical" value obtained from the decay 
A — > 7T + N, i.e., Z^tva/^ 71 " = 0.36. Notice in Table || that the choice of form factor parameterization 
does not have a significant effect on the values of the p and a coupling constants. However, the 
choice of the 7T-/VA vertex does make a difference to k p and gavngaNN ■ It has been shown that A 
pole diagrams constructed using the conventional and Pascalutsa ttNA vertices differ by a contact 



term [13|; in our models it appears that this contact term is partially being mimicked by the p and 



a exchange diagrams. 
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Table 2: Scattering lengths and volumes obtained from the BSE in units of m n , compared 
to results from the SM95 111]] and KH80 [|l^] ttN partial wave analyses. The model I form factor 
parameterization was used. 
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Figure 1: Phase shifts obtained from 
the BSE using the conventional 
ttNA vertex [model I (solid line), 
model II (dot-dashed line)], and the 
Pascalutsa 7riVA vertex [model I 
(dashed line), model II (dotted line)]. 
The data points are from the SM95 
partial wave analysis |TT|. 
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4 Approximations to the Bethe-Salpeter equation 

The BSE has not often been used in meson-exchange models to describe meson-baryon or baryon- 
baryon scattering processes. Due to the relative-energy integration, the BSE is usually regarded as 
being too hard to solve, so it has been much more common to use various approximations to the 
BSE. The X-matrix approach is the simplest method, where a unitary T matrix is found directly 
from the on-shell potential via 



T = V - ilm[G nN ]V. 



(7) 



This is obtained from the BSE when the principle-value parts of all loop diagrams are neglected. In 
Figure 2 we compare the X-matrix approach to the BSE, with the parameters obtained from the 
BSE fits used in the X-matrix approach calculations. For P13 and P31 the if-matrix results agree 
fairly well with the BSE. The results for the remaining partial waves illustrate the importance of 
dressing and multiple scattering: the i^-matrix results deviate significantly from the BSE results. 
Note that for the model II form factors the /f-matrix approach does slightly better than for model 
I, i.e., dressing and multiple scattering are more important in model I than in model II. This is due 
to the large size of the cutoff masses in the model I fit as compared to the cutoff mass in the model 
II fit (see Table §. 

Another approach is to approximate the relative-energy integration in the BSE in some way, 
resulting in a 3-D quasipotential equation. Note that many quasipotential equations depend on the 
choice of fx and fi n [as defined in Eq. (2)] due to the violation of Lorentz-invariance. Here we 
consider the usual choice 
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(8) 



In the Cohen equation [14] it is assumed that the T matrix is independent of the relative-energy, 
and so the relative-energy integration can be performed explicitly over V and G^n, resulting in a 
3-D equation. In Salpeter's instantaneous equation J15[| , it is assumed that the interaction kernel is 
independent of the relative-energy, hence allowing the relative-energy integration to be performed 
explicitly over just G n N- There are an infinite number of possible 3-D equations which are obtained 
from the BSE by replacing G n N by an approximate 2-body propagator which generates the irN 
unitarity cut, but contains a ^-function on the relative-energy. The Blankenbecler-Sugar (BbS) 




Figure 2: Comparison between the 
BSE [model I (solid line), model II 
(dot-dashed line)], and the K-matrix 
approximation [model I (dashed 
line), model II (dotted line)]. The 
conventional ttN A vertex was used. 
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equation [|T(| uses the following propagator: 

G?W(q; P) = 2vri f°° [m f + rf + m N }6^[( m P' + qf - [(/^P' - qf - m% 

J Sth s s le 

(9) 

where s t h = (tun + m^f, and the function f(s', s) is taken as unity. We note that the choice of / 
is in fact arbitrary, provided f(s, s) = 1, hence allowing for the possibility of an infinity of different 
equations. Different choices for /ijy and fi n also result in different equations. The Cooper-Jennings 
(CJ) equation [17] makes use of the propagator 



G&fa P) = 2ni[ m P+i + m N ] d f f{s '' s) 5 M [(VNP' + <1? ~ m%]5^[(^P' - qf - m\\ 

L , s — s — le 



where / is chosen such that G^i, can be rewritten as 



Gffif(q; P) = 2vri ^ [fi N P+ i + m N ] 



(10) 



(11) 
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Using the parameters obtained in the BSE fits to the ttN partial wave analysis, we calculate 
phase shifts using the four different quasipotential equations and compare the results to the BSE. 
In Figure 3 we see that all four of the considered 3-D equations agree reasonably well with the BSE 
in the 531, Pi 3 and P31 partial waves. In the P\\ partial wave, the Cohen and instantaneous equa- 
tions agree well with the BSE results, but the Blankenbecler-Sugar and Cooper-Jennings equations 
generate far too much attraction. For P33 only the Cohen equation gives phase shifts with the same 
shape as the BSE results. The other three equations produce so much attraction that the A(1232) 
resonance has become a bound state, therefore causing the phase shifts at the ttN threshold to be 
180°, rather than 0°. 

The agreement between the Cohen equation and the BSE is not surprising, and has been found 
before in <fxj) scattering [18|. The relative-energy integration over VG^n produces 4-body as well as 
2-body thresholds. Consequently, the T matrix obtained from the Cohen equation contains more 
of the analytic structure generated by the BSE as compared to the other quasipotential equations 
considered here, which only contain the 2-body threshold. 
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Concluding remarks 



'■i ; .. £ Figure 3: Comparison between the 

* < Bethe-Salpeter (solid line), Cohen 
Ps1 (dot-dashed line), Salpeter (long- 

dashed line), Blankenbecler-Sugar 
(short-dashed line), and Cooper- 
Jennings (dotted line) equations. 
The conventional 7riVA vertex was 
/ used, with the model I form factor 

/ f parameterization. 
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In summary, we have presented a relativistic description of 7ryY scattering based on the Bethe- 
Salpeter equation, which we have compared to some approximations schemes, including the K- 
matrix approach and four different 3-D quasipotential equations. In some partial waves large dif- 
ferences were found between the BSE and the other approaches. 

The present model could be extended to higher energies by including resonances into the 
interaction kernel and couplings to inelastic channels, or extended to include the photon so as to 
give a covariant model of pion photoproduction including final-state interactions. Work in both 
these directions is in progress. 
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